Instabilities of micro-phase separated Coulombic systems in constant electric fields. 
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Mixtures of near-symmetric oppositely charged components with strong attractive short range 
interactions exhibit ordered lamellar phases at low temperatures. In the strong segregation limit the 
state of these systems can be described by the location of the interfaces between the components. It 
has previously been shown that these systems are stable against small deformations of the interfaces. 
We examine their stability in the presence of a uniform electric field. When the field is perpendicular 
to the lamellae, the system is unstable against long wavelength deformations for all non-zero values 
of the external field. A field parallel to the lamellae produces deformed but persistent interfaces. 
In a finite thickness system, onset of an external perpendicular field modifies the ground state. 
Flow between the old and new ground states requires the destruction of the original interfaces; this 
destruction proceeds through the instabilities identified in the bulk case. We examine the possibility 
of dynamical stabilization of structures by means of oscillating fields. 



In binary mixtures, the competition between repul- 
sive short-range, and attractive long-range interactions 
leads to micro-phase segregation and pattern formation. 
A basic example is a mixture of immiscible, oppositely 
charged particles, where electrostatic interactions frus- 
trate the full phase segregation of the species. This basic 
scenario has been considered in contexts such as charged 
polymer blends , polyelectrolytes , multi-component 
micelles la| , alloys jj, 0| , and photostimulated semicon- 
ductors |fj . These systems can be studied at a more 
abstract level 0,0 since concrete models often turn out 
to be equivalent to those describing other systems such 
as multiblock copolymers [jj, ^} and reaction-diffusion 
systems 0, i among others. 

Control or modification of self-assembled structures is 
of course of great practical importance. It is natural to 
consider external electric fields as a mean of morphology 
control [l^, as these directly couple to the local properties 
of the system and are directly implementable. Morphol- 
ogy manipulation by electric fields has already proven 
successful in block copolymer systems in which the mi- 
crophases are not charged but have different dielectric 
constants. 0, flil flfij | We investigate in this article the 
effects of coupling the charged binary mixture with an 
external electric field. We describe the ground state of 
the system in the presence of the field; this turns out to 
also be a lamellar state with suitably modified thickness 
at the edges of the system. Our main result is that the 
original lamellar geometry, at the onset of a perpendicu- 
lar external field, is an unstable equilibrium state of the 
system. We use this result and the nature of the pertur- 
bation created by external fields parallel to the lamellae 
to describe some of the main features of the dynamics 
that bridges the original and final equilibrium state of 
the system. 

The model we discuss operates in the strong segre- 



gation limit where concentration fluctuation effects are 
negligible, and where the electrostatic interaction is weak 
compared to the segregation tendencies of the mixture. 
At each point in space there is a well defined major- 
ity component and charge. Between these oppositely 
charged regions there are sharp interfaces and there is 
a well defined surface tension 7 between them. This ap- 
proximation is useful when the characteristic thickness of 
the interface is much smaller than the characteristic size 
of the segregated domains. The segregated regions are 
considered large with respect to the particle size and can 
be treated as continuous media with prescribed charge 
density po. 

We can write an effective free energy for the system 
with contributions from the interfacial surface tension, 
the electrostatic self-energy, and the interaction with the 
external field. We consider only the symmetric case 
where the magnitude of the charge density p is the same 
for each of the species. Dimensional analysis shows that 
the system has characteristic length and energy scales 
given by I — (Po/7) 1 / 3 , and g — jl 2 . All physical vari- 
ables are adimentionalized using these scales. We use a 
charge indicator field p that only takes the values ±1. 
The electric field and potential generated by the charges 
in the system are E; and fa, while E e and <f) e correspond 
to the externally imposed field and its associated poten- 
tial. Using these conventions, the free energy of the sys- 
tem is 
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The region occupied by the charges is labelled A, and the 
interfaces within the system are labelled by T. In the last 
term we introduce a dimensionless pressure field ip that 
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FIG. 1: Panel (a) shows the ground state lamellar structure, 
with layers of width A . White regions are positively charged, 
shaded regions are negative. The scheme shows our choice of 
unit cell and coordinate system and the assignment of labels 
r + , r~ to the interfaces. Panels (b-d) show different normal 
deformation modes; (b) is the charge antisymmetric or wiggle 
mode, (c) is the charge symmetric or corrugation mode, and 
(d) appears at the edge of the Brillouin zone. Panel (b) also 
shows the surface displacement vector u. 



acts as a Lagrangian multiplier enforcing the condition of 
incompressibility. This term does not contribute to the 
net value of the energy but restricts the possible confor- 
mations of the system and its allowed deformations. 

For components with equal or near equal charge den- 
sity, the ground state of the system has a lamellar ge- 
ometry with alternating layers of homogenous charge as 
shown in Fig. In the symmetric case, in zero exter- 
nal field, the net bulk free energy per unit volume / for 
a lamellar pattern of layers of width A is 
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The optimal value of the thickness is then Ao = (3/tt) 1 / 3 . 
The internal electric field is not zero, and generates a 
non-zero force within the bulk of the system. This force 
is compensated by the pressure gradient V'0 and, in fact, 
this condition determines the values of the pressure field. 

Analysis of the stability of the lamellar ground state of 
this system in the absence of external fields has been ana- 
lyzed by Muratov Q • Below we sketch, in notation useful 
for our problem, a derivation of the stability properties 
of the system, ending the main result of Eq|3| Later, 
we use the results in the context of external fields. The 
fluctuations of the system can be described by specifying 
the deformations of the interfaces. We use a coordinate 
system x = (r, z), where the vector r runs parallel to 
the lamellae and the coordinate z is transversal to them. 
The unit cell of the one-dimensional lattice formed by the 
lamellae can be taken to be symmetric with respect to its 
mid-plane. We assign z = to the midplane so that the 
cell extends up to z — ±Aq while the interfaces, labelled 



FIG. 2: Ground states of a finite thickness system. On the 
left, in the absence of an external filed, the system is symmet- 
ric. The thick lines indicate the boundaries of the system. If 
the unit cell width 2Ao is commensurate with the thickness, 
the outmost layers have equal charge and thickness Ao/2. As 
shown in the right frame, in the presence of an external field 
the ground state is obtained by shifting the location of the 
interfaces. The new charge distribution at the boundary com- 
pensates the external field in the bulk. 



r , are located at z± = ±Ao/2. The collective deforma- 
tions of the system prescribe a transversal deformation 
A u n,± = u n ,±z for each of the interfaces r„ , where n is 
an index for the cell. The normal modes of the system 
have the form u n ,± = A±cos(2pn\o) cos(q • r), where the 
wave-vector in the transversal direction q is unrestricted 
while the longitudinal wavevector p takes values within 
the first Brillouin zone \p\ < tt/(2\q). There are two 
normal modes for each wave vector (p, q) , each with a 
different ratio between the deformation amplitudes A±. 

A small deformation of the interfaces has the effect of 
increasing the interfacial area, and producing changes in 
the charge density localized at the interface. The ex- 
cess charge density interacts with the background charge 
and with itself. The excess surface tension energy den- 
sity associated with a small deformation has the form 
(l/2)(V||u(z)) 2 where V|| is the gradient along the inter- 
face. The excess charge density can be formally written 
as a series expansion in powers of the deformation ampli- 
tude. The terms of these series, pW, and the associated 
series for the the potential and field are labelled by i, the 
order of the term in the expansion. The first terms of the 
excess charge density at a given interface have the form 
P = P^ T 2u(v)8(z — z±) ± u(r) 2 S'(z — z±) where 8 and 
8' are the delta function and its derivative. The poten- 
tial and electric field can be solved order by order in the 
series. The higher order terms of the potential satisfy 
the Laplace equation within each homogenous domain 
and have boundary conditions determined by the excess 
charge density at the interfaces. The electrostatic excess 
energy density is {l/2)4>^ + p^(/)^ \ Upon integra- 
tion of the delta functions localized at the interfaces, the 
excess free energy is the sum, over all interfaces, of the 
terms 



AF„,± = / dr 



(V|,u) 2 + . 
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We choose the orientation of the interface so that its nor- 
mal vector n points outward from the positive charge re- 
gions. The electric field at the interfaces in the ground 



3 



state is E^ -* = 27rAon, so that the last term in the excess 
energy always has a destabilizing effect. 

Carrying out the calculation sketched above, the en- 
ergies of the normal modes, as well as the relation be- 
tween the deformation amplitudes of the surfaces in the 
unit cell can be obtained as a function of the wavevector 
(p, q) . An equivalent calculation has been carried out by 
MuratovQ, including the case of asymmetric mixtures, 
and in the presence of screening. The basic result is that 
all deformations have a non-negative excess energy and 
therefore the system is stable against small deformations. 
As shown below, in the presence of constant external 
fields some of the modes become unstable. The modes 
that become unstable are deformations of the lowest en- 
ergy modes at zero field that are located at the center, 
at p = 0, and edges p = ±7r/(2Ao) of the first Brillouin 
zone. The explicit form and energies of these modes at 
zero field are as follow. For p — 0, all unit cells have 
the same displacements, and within each unit cell the 
normal modes are the charge symmetric A± = ±A s /v2, 
and charge antisymmetric modes A± = A a j\[2. The en- 
ergy per unit volume of each of these modes, for each 
transverse wavevector q with magnitude q, are (l/2)eA 2 
with 
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These results exhibit the softness of the antisymmetric 
mode, e s ~ q 4 and the divergent hardness of the symmet- 
ric mode e s ~ q~ 2 , for small q. In the edge normal modes 
with p = ±7r/(2Ao) only one of the unit cell surfaces is 
deformed. The two modes correspond to the choice of 
deformed surface and clearly both modes have the same 
energy density (l/2)e e A 2 with: 
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The lowest energy modes in this group appear at gAo = 
1.6. 

Finite size effects and the action of external fields are 
intimately related in this system. Before addressing the 
external field, it is useful to consider the properties of 
the ground state of the finite system. First, we point 
out that the ground state of a slab of material of width 
commensurate to 2Ao is identical to the bulk solution, 
with lamellae parallel to the limiting surfaces, terminat- 
ing at both ends with layers of the same charge of size 
Ao/2. Consideration of the surface terms generated by in- 
tegration by parts during variation of the charge density 
leads in principle to boundary conditions for the system. 
In this case, however, these turn out to be equivalent to 
the bulk equilibrium conditions, and the global minimum 



must be determined from evaluation of the bulk energy 
restricted by the incompressibility condition but without 
extra boundary requirements. 

For non-commensurate thicknesses the energy mini- 
mum is obtained when the first two layers at both ends 
are still arranged symmetrically but have thickness differ- 
ent from Ao- If the end layers do not have equal charge, a 
net electric field is created within the system. The ground 
state of the finite system in the presence of a constant 
perpendicular field is then simply a lamellar structure 
with unbalanced end layers that create an excess inter- 
nal field that exactly compensates the external field, as 
sketched in Fig. 

We consider now the effect of the sudden onset of an 
external uniform electric field on the ground state at zero 
field. We first note that the original lamellar structure is 
still an equilibrium state even after the onset of an ex- 
ternal field perpendicular to the lamellae. Indeed, due to 
the incompressibility of the system, rigid translations of 
the interfaces between components are not possible and 
all other deformations of the interface lead, to first order, 
to net zero changes in the energy of the system. In more 
detail, the equilibrium of the system is achieved again by 
letting the internal pressure gradient to adjust to com- 
pensate the net electric field. This effect can always be 
achieved when the external field has null curl, as is the 
case of an uniform field. As we show below, however, this 
equilibrium state is unstable. 

The external field E e ± perpendicular to the lamellae 
couples to the second order term of the charge distribu- 
tion induced by deformations. The excess energy is given 
by the same expression as the one obtained for the inter- 
action of the fluctuations with the internally generated 
field in Eq. J3J). For a given surface, this energy is 



AF = - 



u 2 E 



n. 
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This interaction has a destabilizing effect only in the in- 
terfaces whose normal is parallel to the external field. 
When the deformations are expressed in terms of zero- 
field normal modes, the field couples to a mixture of the 
p = symmetric and antisymmetric modes. The excess 
energy density is A/ = E e± (A_{q) 2 - A + (q) 2 )/(4A ) = 
— A s A a E e ±/ X . The new p = normal modes in the 
presence of the field are obtained from diagonalization 
of the new energy quadratic form. These new modes 
are smooth modifications of the zero-field modes, and 
will be labelled by the corresponding zero-field labels 
a or s. While it is clear that non-zero average excess 
energies arise for p = modes, the coupling to the 
square of the amplitude also gives non-vanishing contri- 
butions for the modes at the edge of the Brillouin zone at 
p = ±7r/(2Ao). The external field shifts the edge modes 
by A/ = ±A 2 E e± /(2\ ). 

There are two sets of unstable modes in the system 
in the presence of an external perpendicular field. The 
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FIG. 3: Plots of the coefficients of the energy quadratic form 
as a function of the transversal wave vector magnitude q. Re- 
sults for the p — charge antisymmetric modes are shown on 
the left, and Brillouin zone edge modes with p = 7r/(2Ao) on 
the right. The labels correspond in both cases to the magni- 
tude of the external transversal electric field (a) E± — 0, (b) 
E± = 0.5,(c) E± = 1.0, and (d) E± = 1.5. Unstable modes 
are those with e < 0. 



modes obtained from the deformation of the antisym- 
metric charge mode, with p — 0, become unstable for 
any non-zero value of the external field. These unsta- 
ble modes appear within a wavevector disk, centered at 
q = 0, of radius q r ~ E e . The most unstable wavevector 
has q m ax ~ E e and the energy coefficient of this mode 
scales as e ~ —E*. Therefore at very low field strengths 
the instability is very weak while, at stronger fields, the 
dynamics are dominated by the exponential growth of fi- 
nite wavelength modes. A band of edge modes around 
the softest edge mode, located at p — it/ (2Ao), q se = 1.6, 
also becomes unstable at the finite field value E ec = 0.9. 
Above this threshold value the band of unstable modes 
has width Aq ~ {E—Eec) 1 ' 2 . The most unstable mode is 
always q se , with energy coefficient e e (q ec ) — E/Xq. Plots 
of the energy coefficients for these modes appear in Fig. 

El 

We consider now the effect of a field Em parallel to 
the lamellae. The ultimate end effect of this field will 
be to orient the lamellae in a direction perpendicular to 
the field. There are, however, dynamical steady states 
with orientations parallel to the field. To show this, 
we decompose the external potential into its Fourier 
components </> q . These components couple linearly to 
the symmetric deformation modes of same wavevector 
and their excess energy density in this case is A/ = 
(l/2)e s (q)A s (q) 2 + V20 q A s (q)/A o . Therefore, start- 
ing from a parallel orientation, the system flows to a 
metastable state with symmetric deformations of ampli- 
tude .A s (q) = V2(t)q/(\o£s(q))- In a system where the 
external field is locally constant but periodic, with pe- 
riodicity L in the direction of the lamellae (see Fig. |3J, 
the potential has non-vanishing Fourier components of 
the form </> q ~ E\\q~ 2 . We pointed out above a simi- 
lar behavior for the self-energy of the long wavelength 
symmetric fluctuations £ s ~ q 2 - The deformation am- 
plitudes are then near uniform for small wavectors, and 
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FIG. 4: Effect of uniform fields parallel to the lamellae, (a) 
Periodically alternating uniform fields cause a concentration 
of charge at the ends of the uniform region, (b) In an infinite 
system, steady currents are possible; the arrows indicate flow 
direction, (c) In finite systems the field induces smooth flows, 
(d) Reversal of the field simply reverses the flow but does not 
destabilize it. 

very small for large ones. The net result for large peri- 
odicities L is an accumulation of charge at the edges of 
each homogeneous region as in Fig. QJi. These deforma- 
tions require net charge exchanges between the bound- 
aries of the homogenous field regions. In the limit of 
infinite size, this charge redistribution emerges as steady 
exchange currents as sketched in Fig. 0Jj. 

The ground states corresponding to different external 
fields cannot be connected by charge redistributions that 
preserve the integrity of the layers since there are no 
smooth deformations that change the charge of the end 
layers. Therefore, the flow between ground states always 
occurs by means of the destruction of the original struc- 
ture by the growth of instabilities. The flow between 
the ground states for different external fields requires net 
transport of charge which, as noted above, occurs natu- 
rally in lamellar orientations parallel to the external field. 
After the destruction of the layers parallel to the bound- 
ary, charge transfer might occur through the formation 
of layers perpendicular to the boundary but parallel to 
the field, until the charge unbalance necessary to shield 
the bulk from the external field is achieved. Afterwards, 
the intermediate layers may become unstable and flow 
into the final stable conformation with layers parallel to 
the boundary. 

While the presence of instabilities precludes the 
smooth flow between ground states, i.e. adiabatic con- 
trol of the system by tuning of external fields, there 
are interesting possibilities for the creation of structures 
through slowly varying uniform fields. In particular, we 
note that smooth deformations of lamellae under paral- 
lel fields carry out a similar process as a growing unsta- 
ble mode in a perpendicular field; namely, both transfer 
charge along the field direction. That is, from an en- 
ergetic point of view, the preferred intermediate states 
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between the ground states of a slab with perpendicular 
fields are those that connect same-charge regions across 
the slab. When driven by an oscillating uniform field, the 
system may oscillate along a set of states of easy charge 
transfer, as sketched in Fig. BJc-d), within a suitable 
range of frequencies determined by the specific dynamics 
of the system. In other words, oscillating fields may dy- 
namically stabilize lamellar orientations parallel to them. 

In conclusion, we have shown how to use the static 
stability analysis of the system as a guide to describe 
possible dynamic scenarios for the flow of the system to- 
wards equilibrium. The key result is that lamellar struc- 
tures are unstable equilibrium states in the presence of 
external perpendicular electric fields, while fields parallel 
to the lamellae give rise to natural charge transferring 
currents. 
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